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Method  of  Approach 


Sediment  transportation,  .like  many  another  engineering  prob- 
lem, has  long  been  studied  in  a  purely  empirical  manner  in  the 
attempt  to  develop  useful  working  rules  from  a\  limited  supply  of 
measured  data.    Since  formulas  for  predicting  the  sediment  load  of 
a  stream  are  generally  based  upon  past  records-  for  similar  types  of 
flow,  they  can  be  depended  upon  to  yield  satisfactory  results  only 
within  this  range  of  flow  conditions.    Nevertheless,  any  empirical 
formula  is  almost  certain  to  be  extrapolated  f 'ar  beyond  its  limit 
of  validity — a  procedure  leading  sooner  or  later  to  serious  error, 
for  it  is  practically  impossible  to  obtain  and!  evaluate  enough 
field  measurements  to  establish  empirically  a  ttruly  general  trans- 
port law. 

Purely  rational  methods,  on  the  other  basnd,  are  even  less 
likely  to  provide  an  ultimate  solution,  owing  to  the  many  simplifi- 
cations which  are  necessary  before  a  rigorous  -mathematical  treat- 
ment becomes  possible.    In  other  words,  neither  mathematical  tools 
nor  physical  understanding  of  their  use  can  be  considered  suffi- 
ciently far  advanced  to  cope  with  so  intricate  a  problem  at  the 
present  time. 

In  recent  years,  however,  considerable  progress  in  hydraulic 
research  has  resulted  from  methods  that  are  a  "healthy"  combination 
of  the  empirical  and  the  rational,  the  process;  of  physical  analysis 
being  closely  correlated  with  experimental  observation.    That  is, 
after  eliminating  those  factors  which  observation  has  shown  to  be 
of  minor  importance^  the  essential  features  of  the  problem  are  an- 
alyzed in  as  rational  a  manner  as  possible.    All  conclusions  are 
then  subjected  to  a  thorough  experimental  check,  the  number  of  data 
necessary  to  the  purely  empirical  approach  hawing  been  greatly 
reduced  through  systematic  organization  of  the  variables. 
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The  practicability  of  this  method  is  evidenced  by  the  prog- 
ress which  it  has  made  possible  in  the  study  of  fluid  turbulence. 
Although  turbulent  flow  is  still  not  fully  understood,  means  have 
been  found  of  estimating  the  velocity  distribution  and  resistance 
of  flow  along  smooth  and  rough  boundaries  in  a  manner  which  has 
possibilities  of  very  broad  application.    Not  only  should  such 
methods  prove  equally  useful  in  the  study  of  sediment  transporta- 
tion, but  existing  knowledge  of  turbulent  flow  should  have  direct 
bee-ring  upon  this  problem.    Indeed,  the  analysis  of  sediment  sus- 
pension purely  as  a  phenomenon  of  fluid  turbulence  has  already 
attracted  considerable  attention. 

For  a  number  of  reo.sons,  the  movement  of  bed  load  and  the 
movement  of  suspended  load  have  always  been  treated  as  two  unre- 
lated phases  of  the  transport  problem.    In  many  localities  the 
sediment  load  of  streams  is  for  the  most  part  feoo  coarse  to  be 
carried  in  suspension,  and  the  problem  remains-  one  of  bed  load 
alone.    In  other  instances  the  task  of  the  engineer  has  been  to 
keep  a  watercourse  open  to  traffic,  and  the  regulation  of  the 
coarser  material  moving  near  the  bed  has  been  Bais  sole  concern. 
These  circumstances,  in  turn,  have  focused  attention  upon  the 
hydraulic  river  model,  in  which  only  the  bed  l:cad  can  be  properly 
simulated.    As  a  matter  of  fact,  not  until  recent  years  has  sus- 
pended material  begun  to  receive  its  share  of  attention  because  of 
the  growing  importance  of  silt  control  in  reservoirs  and  irrigation 
projects . 

While  a  fair  degree  of  success  has  accompanied  the  study  of 
bed  load  by  purely  empirical  methods,  these  Yimve  not  seemed  to  lend 
themselves  so  readily  to  the  more  advanced  staige  of  sediment  sus- 
pension.   On  the  other  hand,  though  the  suspension  of  sediment  by 
turbulence  is  apparently  more  susceptible  to  r-rational  analysis,  the 
treatment  of  suspended  load  as  a  completely  independent  problem  has 
proved  to  be  Impossible,  for  the  turbulence  of  the  flow  and  the 
amount  of  material  carried  into  suspension  are  both  governed  by 
conditions  in  the  bed  region.    It  is  evident,   therefore,  that  an 
analysis  of  the  sediment  problem  as  a  whole  wiill  first  become  pos- 
sible when  bed  load  and  suspended  load  can  be  expressed  as  func- 
tions of  the  same  flow  parameters. 

Turbulence  and  Bed  Load. 

Although  flow  boundaries  may  be  loosely  .classed  as  either 
rough  or  smooth,  boundary  roughness  is  actually  a  relative  matter. 
It  has  been  found  expedient  in  the  study  of  fl.ow  through  pipes  (9)^ 


'Numbers  in  parentheses  refer  to  Literature  Cited,  pp.  24-25. 
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to  define  relet ive  roughness  as  the  ratio  of  a  characteristic 
roughness  height  to  the  radius  or  diameter  of  the  pipe;  even  then, 
however,  a  pipe  of  a  given  relative  roughness  may  be  effectively 
smooth  at  one  Reynolds  number,  and  rough  at  another.    vVhile  the 
relative  roughness  of  the  pipe  and  the  Reynolds  number  of  the  flow 
are  the  two  parameters  which  determine  the  type  of  flow  to  be  ex- 
pected, considerable  light  is  thrown  on  the  problem  by  introducing 
another  type  of  relative  roughness — the  absolute  roughness  of  the 
boundary  in  its  ratio  to  the  thickness  of  the  zone  of  laminar  flow 
along  the  pipe  wall. 

Within  the  thin  laminar  layer  at  a  very  smooth  boundary 
(shown  at  greatly  exaggerated  scale  in  fig.  1),  the  velocity  in- 
crease? so  rapidly  from  aero  to  that  of  the  surrounding  turbulent 
flow  that  it  may  bo  assumed  to  vary  with  the  first  power  of  the 
distance  from  the  boundary  surface.    Thus,  in  terms  of  the  inten- 
sity of  boundary  shear  and  the  dynamic  viscosity  of  the  flui  d3 

or,  introducing  the  density  and  kinematic  viscosity,  and  rewriting 
in  dimensionless  form 

v  _  ffo  y  (2) 


The  factor/^//5     has  the  dimension  of  a  velocity,  and  is  commonly 
known  as  the  "friction  velocity."    It  is  proportional  to  the  mean 
velocity  of  flow  and  the  square  root  of  the  resistance  coefficient, 


VZ/p  =  VVf/s 


In  the  zone  of  fully  developed  turbulence  in  very  smooth 
pipes,  velocity  measurements  follow  the  analytical  relationship 


(3) 


If  the  thickness  of  the  boundary  layer  is  arbitrarily  defined  as 
the  intersection  of  the  velocity  curves  of  equations  2  and  3 
(fig.  1),  it  will  be  found  that 


3?or  table  of  symbols  and  their  significance  refer  to  Tabulation 
of  Symbols  on  pp.  22-23- 
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whereupon  equation  3  may  be  rewritten  in  terms  of  boundary- layer 
thickness : 

=    //.6  +  5'.7S  loo  JL.  (5) 


00 

d 


Nikuradse  (6),  vrho  determined  the  constants  of  equation  35 
extended  his  research  to  pipes  which  had  been  artificially  rough- 
ened by  cementing  sand  grains  on  the  inner  surface.    Using  the 
grain  diameter,  d  ,  as  a  measure  of  absolute  roughness,  k  ,  the 
velocity  distribution  in  such  pipes  was  found  to  be  expressible  in 
the  form 

v    -      u  /nzr      1  1 

(6) 


in  which  the  function  ^y<F^b  -—J       will  be  seen  to  be  equiva- 
lent to  (p'fc/ f  fj  .    In  other  words, 

 5.75  iog  4  ,  )  (V) 

Yt/p  *   '  '  s  J 

A  plot  of  measured  data  in  the:  form     v  l/ 71  / p     ~  5 .75  log  y 

versus   a/ /$    is  shown  in  figure  2j  the  ratio  of  absolute  roughness 
to  boundary-layer  thickness  is  seen  to  be  of  decided  importance. 

For  very  low  values  of  or  /  £    ,  the  wall  irregularities  are 
fully  enclosed  by  the  laminar  layer  (fig.  3a),  and  the  flow  is  then 
effectively  the  same  as  that  through  an  absolutely  smooth  pipe. 
Equation  5,  for  smooth  pipes,  should  describe  such  motion  when  re- 
written in  terms  of  roughness  parameters.    Subtracting  5 • 75  log  Q /c/ 
from  both  sides  of  this  expression,  it  will  be  seen  that 

-  5-75  log-^    =  5.75  log  -JSi  +  11.6  (8) 
c/  S 


ix/p 

which  is  the  equation  of  the  sloping  line  at  the  left  of  figure  2, 
approached  asymptotically  by  the  measured  data  as  d / $  becomes 
very  small.    At  moderate  values  of  c/ / S      (fig.  3b),  the  boundary 
irregularities  begin  to  influence  the  flow  outside  the  laminar 
film,  the  boundary  then  being  no  longer  effectively  smooth.  With 
increasing  values  of  c// J    ,  the  turbulence  induced  by  the  wall 
roughness  becomes  greater,  until  finally  the  flow  at  the  boundary 
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is  fully  turbulent  (fig.  3c).    Under  such  circumstances,  the  func- 
tion y(\/Z^Jp  o//^~  <p>  '('cJ/'s)     °^  equation  7  becomes  constant, 

since  the  velocity  distribution  is  no  longer  dependent  upon  viscos- 
ity. Thus,  as  will  be  seen  from  figure  2,  for  very  large  values  of 
a/ /g~  the-  experimental  points  approach  the  horizontal  asymptote 

(^JvJS  .  '  .  ; ,  *  ~~ .  ..  -  -.i  .  _.        ......  |j 

II 

1  v      -  5.75  log         =  8.48  !  (9) 

of  I 
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Because  of  the  interdependence  of  velocity  distribution  and 
boundary  resistance,  it  may  also  be  shown  (7)  that 


-  5.75  log  ~  =   ^jg  -  5.75  log  +  3.75         '  (10) 

tZ/p  ■     Yf  d 

the  resistance  factor  for  very  smooth  boundaries'  varying  in  the 
form 

Vf  ' 

while  for  fully  effective  boundary  roughness, 

/  .  K 

rzrzz.    =    1.74  +  2  log  -12—  (12) 

Figure  2  is  thus  seen  to  include  two  series  of  experimental  data, 
one  determined  on  the  basis  of  velocity  measurements,  the  other 
from  measurements  of  resistance.    The  close  agreement  between  the 
two  series  is  good  evidence  of  the  validity  of  both  relationships, 
the  experimental  points  following  a  single  functional  trend  between 
the  limits  of  smooth  and  rough  boundary  flow. 

Since  the  artificial  roughness  used  in  these  experiments  had 
essentially  the  same  texture  as  that  of  a  leveled  sand  bed  in  a 
laboratory  flume,  the  foregoing  relationships  should  be  directly 
applicable  to  the  study  of  the  beginning  of  bed-load  movement. 
This  method  of  attack  was  adopted  by  Shields  (12),  who  reasoned 
that  the  initial  movement  of  material  of  any  given  particle  size 
would  require  the  existence  cf  a  certain  critical  velocity,  vc  ,  at 

a  distance  above  the  bed  proportional  to  the  sediment  diameter, 

Yc     =c(  c/    .    Thus,  from  equation  7, 


Vf  [ 5-75  iog  «-  *  ?(f)\  -  '0  ?•(«. ,  f  )  (l3) 
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The  force  exerted  by  the  flow  under  these  conditions  is  expressible 
in  the  form 

in  which  .1>  is  a  coefficient  of  resistance  to  flow,  depending  in 
magnitude  upon  the  sediment  characteristics  end  a  Reynolds  number 
for  local  flow  conditions j  that  is, 

t  rj«*  >  --f)  (15) 

Introducing    vc    from  equation  13,  X 


The  resistance  to  movement ,  on  the  other  hand,  will  vary  with 
the  effective  weight  of  the  immersed  sediment  and  with  the  form  and 
arrangement  of  the  grains : 

r  =  s(^-^V3  (is) 

At  the  beginning  of  motion,  it  is  evident  that    R    must  equal    F  , 
v/hence. 


in  which    Tc    has  its  customary  significance  of  critical  tractive 
force,  i.  e.,  the  intensity  of  boundary  shear  at  the  beginning  of 
motion.    Since  cC,  f  oci     and     cc^      all  depend  upon  the  character- 
istics of  the  bed  material,  each  should  be  expressible  in  terms  of 
parameters  representing,  statistically,  particle  shape,  grading, 
and  compaction.    Under  such  conditions,  d    properly  refers  to  a 
geometric  mean  size  of  grain.    On  the  other  hand,  if  the  materials 
to  be  studied  are  essentially  uniform,  in  size,  shape,  and  compac- 
tion, equation  19  becomes  simply 
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Shields  conducted  a  series  of  experiment- s  on  the  initial 
movement  of  uniform  materials  having  a  wide  range  in  specific 
weight,  plotting  results  according  to  equation!  20  as  shown  in  fig- 
ure 4«    At  once  apparent  is  the  great  similarity  of  the  functional 
trend  to  that  of  measurements  on  rough  pipes  In  figure  2.    In  other 
words,  the  initial  movement  of  bed  particles  is  directly  influenced 
by  the  conditions  of  motion  in  the  neighborhood  of  the  bed.  The 
4-5°  line  at  the  left,  representing  undisturbed  laminar  flow  at  the 
bed,  has  been  fatted  by  eye,  and  follows  the  equation 


 2L_  =  o.oi 


(21) 


Similarly,  the  horizontal  line  at  the  right  representing  fully  de- 
veloped turbulence  at  the  bed  is  given  by  the  equation 

— =    0.06  (22) 

It  is  obvious  that  initial  movement  may  occur-  within  the  laminar 
boundary  layer,  but  that  the  influence  of  viscosity  upon  such  move- 
ment steadily  decreases  as  the  boundary  layer  is  broken  up  by 
grains  of  increasing  relative  size. 

Of  definite  academic  importance  is  this  study  by  Shields,  in 
that  it  shows  for  the  first  time  the  effect  of  viscous  action  on 
the  beginning  of  movement.    Its  direct  application  is  limited,  how- 
ever, to  the  choice  of  materials  for  model  steadies,  since  beds  in 
natural  streams  display  neither  the  initial  regularity  nor  the 
grain  uniformity  of  the  laboratory  study.    Of  perhaps  even  greater 
importance,  therefore,  is  the  fact  that  Shields  found  the  form  of 
initial  bed  irregularities  to  depend  directly  upon  the  magnitude  of 
c//<f    .    At  lovr  values  of  this  ratio,  for  ins.tanco,  short,  deep 
ripples  developed.    As  the  influence  of  viscosity  decreased,  the 
ripples  tended  to  become  shallower  and  longer,,  approaching  a  limit 
when  the  boundary  layer  was  finally  broken  up'^  completely.  Evident- 
ly, the  use  of  coarse,  low-density  materials  such  as  pumice,  amber, 
and  lignite  to  simulate  bed  sediment  in  river*  models  tends  to  de- 
crease the  influence  of  viscosity  upon  the  bed  development. 

As    T    increased  beyond  the  critical  value    Tc  (following 
the  direction  of  the  arrows  in  fig.  4,  since     T      / /$z  ),  the  bed 

irregularities  were  found  to  increase  in  depth,  but  at  the  same 
time  the  rate  of  bed  movement  became  greater-    Shields  reasoned 
tha.t  the  sediment  transport  should  be  expressible  in  terms  of  dis- 
charge, slope,  and  the  bed  parameter  (T-  7^ }}/{  /f~  </*]  o/  , 
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once  the  motion  at  the  bed  became  fully  turbulent,  namely,  inde- 
pendent of  the  fluid  viscosity.    An  analysis  of  measured  data  (fig. 
5)  proved  that  the  transport  function  could  be  written  in  the 
dimensionless  form 

JL  -   10  s  l~.r<-~.  "  (23) 

the  factor    G    representing  the  rate  of  sediment  transport  by 
weight  (measured  under  water),  and   Q    the  rate  of  water  discharge, 
also  by  weight  (this  relationship  may,  of  course,  be  written  in  a 
number  of  equivalent  forms  through  introduction  of  the  equilibrium 

equation  for  boundary  shear,    T  =  J^DS  ,  and  the  Chesy  equation, 
V  =0*\/£>S  .)    At  low  values  of  (T-Tc)/(  <ts  -  f  )  of }  the  movement 

consisted  of  rolling  or  sliding  along  the  bed,  but  at  more  advanced 
stages  the  sediment  grains  were  carried  over  the  bed  irregularities 
in  the  process  known  as  saltation.    Movement  finally  became  so 
pronounced  that  the  crests  of  the  irregularities  were  v/ashed  away — 

at  approximately  (T-Tc)/(.-^-  <F)d  =  O.2.       .    Few  experimental 
data  are  shown  beyond  this  zone,  for  the  turbulence  of  the  flow, 
once  saltation  developed  to  a  sufficient  degree,  would  carry  the 
saltating  material  into  suspension  (see  fig. 

Equation  23  must  be  regarded  as  essentially  empirical,  in 
particular  since  it  was  based  upon  measurements  with  material  hav- 
ing less  than  a  two-fold  variation  in  size,  and  restricted  to  con- 
ditions under  which  the  material  moved  only  as  bed  load.  Neverthe- 
less, were  the  relationship  sufficiently  well  founded  to  warrant 
extrapolation,  it  could  also  be  used  to  predict  the  total  amount  of 
material  carried  by  the  flow  in  more  advanced  stages  of  suspension. 
Although  the  actual  function  may  not  be  linear  oven  in  the  region 
shown  in  figure  5,  there  is  no  reason  to  expect  an  abrupt  discon- 
tinuity at  any  point.    Nevertheless,  the  possibilities  of  develop- 
ing such  a  general  relationship  do  not  seem  to  have  been  investi- 
gated.   Experimental  measurements  invariably  cease  once  an 
appreciable  amount  of  sediment  goes  into  suspension,  due,  in  part 
to  the  difficulty  of  trapping  all  material  in  the  usual  bed-load 
flume,  and  in  part  to  the  custom  of  regarding  sediment  suspension 
as  a  phenomenon  of  totally  different  category. 

Turbulence  and  Suspended  Load 

Owing  to  the  belief  that  flow  conditions  in  the  bed  region 
differ  from  those  at  some  distance  above  the  bed,  previous  studies 
of  the  mechanics  of  suspension  have,  in  turn,  eliminated  the  bed 
zone  from  consideration.    Thus  limiting  the  analysis  to  regions  of 
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fully  developed  turbulence,  it  has  repeatedly  been  shown  that  the 
sediment  concentration  at  any  level  (fig.  6)  shculd  vary  according 
to  the  function 


Herein    c    represents  the  concentration  of  the  suspension  in  per- 
cent by  weight,  w    is  the  terminal  velocity  of  settling  of  the 
sediment  grains  in  stagnant  water,  and  <f  is  a  kinematic  turbulence 
coefficient,  corresponding  roughly. to  the  kinematic  viscosity — that 
is,  while  V  is  a  measure  of  the  intensity  of  molecular  mixing,  £ 
refers  to  a  molar  mixing  process  in  which  finite  masses  of  fluid 
participate  in  eddy  motion.    This  process  of  intermixing  between 
neighboring  regions  of  flow  tends  to  make  the  sediment  distribution 
more  uniform  by  carrying  material  in  the  direction  of  decreasing 
concentration — upwards,  in  the  normal  case — at  a  rate  proportional 
to  the  intensity  of  mixing  and  to  the  rate  of  change  of  concentra- 
tion at  that  level.    Because  of  the  weight  of  the  sediment,  on  the 
other  hand,  material  tends  to  move  downwards  at  a  rate  proportional 
to  the  concentration  at  that  level  and  to  the  velocity  of  settling. 
The  state  of  equilibrium  described  by  equation  24  presumes  the 
effects  of  settling  and  mixing  to  be  in  balance  for  all  values  of 
y    at  which  the  flow  is  turbulent, 

Since  the  mixing  coefficient  is  diniensionally  equivalent  to 
the  product  of  a  velocity  and  a  length,  the  velocity  may  be  con- 
sidered proportional  to    v'y  ,  the  mean  absolute  velocity  of  the 

eddy  currents  in  the  vertical  direction,  and  the  length,  ^  ,  pro- 
portional to  either  the  average  size  of  the  eddies  or  the  vertical 
distance  the  fluid  masses  are  carried  before  losing  their  identity 
in  the  surrounding  fluid.  While  the  fall  velocity,  w  ,  refers  to  a 
specific  size  of  sediment,  <f  must  be  regarded  purely  as  a  statis- 
tical mean  product,  both    v'y    and  Jc    varying  with  tine  in  accord 

with  the  normal  error  law.  It  is  pointless,  therefore,  to  presume 
that    w    m:ust  be  smaller  than  the  mean  absolute  value  of    v'  be- 

fere  suspension  may  exist,  for  there  is  no  limit  to  the  magnitude 
of  instantaneous  fluctuations.    However,  equation  24  does  require 
that  the  mixing  of  the  sediment  carried  by  the  eddies  be  statisti- 
cally the  same  as,  or  at  least  proportional  to,  that  of  the  fluid 
itself.    The  mixing  coefficient  should  thus  be  written  as 


c/c 


=  ~CLV 


(24) 


(25) 


in  which  @    is  a  bulk  coefficient  of  proportionality. 
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The  general  validity  of  equation  24  was  verified  within  the 
past  year  by  means  of  an  agitating  mechanism  producing  a  constant 
intensity  of  mixing  throughout  a  tan]-:  of  water  (11)-    Under  such 
circumstances   <f    becomes  independent  of    y  ,  whereupon  equation 
24  may  be  integrated  to  yield  the  concentration  at  any  level  in  its 
ratio  to  that  at  some  arbitrary  level    a    (see  fig.  6);  thus  - 


Four  different  sizes  of  uniform  quartz  sand  were  tested,  permitting 
an  approximately  30-fold  variation  in  fall  velocity.  ■  The  magnitude 
of  £  varied  directly  with  the  frequency  of  agitation,  which  could 
be  controlled  mechanically  over  an  8-fold  range.    Curves  of  sedi- 
ment distribution  could  be  determined  from  pipette  samples  taken 
over  a  vertical  distance  of  1  1/2  feet.    Not  only  were  the  various 
sizes  investigated  individually,  but  a  mixture  of  all  four  grades 
was  tested  at  various  values  of   £  ,  yielding  t.he  composite  dimcn- 
sionless  plot  shewn  in  figure  7.    Aside  from  secondary  deviations 
(due  to  variations  in  both  Q  and    w    for  the  coarse  material  in 
the  presence  of  fine),  the  points  follow  the  exponential  curve  of 
equation  26  in  a  sr.tisfactory  manner. 

Application  of  equation  24  to  the  suspension  of  sediment  in 
f lowing  water  is  a  more  difficult  matter,  since  <f  then  varies  as  a 
function  of    y  .    For  want  of  knowledge  as  to  its  actual  variation 
in  an  open  channel,  one  is  compelled  to  turn  again  to  studies  of 
flow  in  pipes.    In  deriving  expressions  for  velocity  distribution 
and  resistance  (refer  to  fig.  8),  Prandtl  assumed  that  the  mean 
absolute  component  of  fluctuation  in  the  direction  of  flow,  v'y  , 

would  be  proportional  to  the  product  of  the  mixing  length  and  the 
mean  velocity  gradient  ( V^'  ^~  A  dv /c/y    )}  and  that  the  components 
of  fluctuation  in  the  two  directions  would  be  proportional  to  each 
other  (       ^  v         ) .    Since,  as  shown  by  Reynolds,  the  intensity 
of  shear  varies  with  the  mean  product  of  the  fluctuations 

(  2"  -  ~  f5 .  K  vy         ),  and  since  intensity  of  shear  is  distributed 
linearly  along  the  pipe  radius  according  to' the  expression 

T  =    TQ  (  /  ~y/V  ),  it  will  be  seen  that  the  following  equations 
for    v'y  ,  J?  3  and  <f  (assuming  £    and    v'y    to  absorb  the  pro- 
portionality factors)  should  then  hold: 


c 


=  e 


(26) 


(27) 
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The  variation  of    v'y    appears  to  be  independent  of  the  ve- 
locity gradient ,  and  is  plotted  in  dimensionless  form  as  curve  |  A 
of  figure  9.    This  curve  is  seen  to  extend  from  zero  at  the  axis  to 
a  maximum  at  the  wall.    However,  Prandtl  did  not  intend  equations 
27,  28,  and  29  to  apply  at  either  of  these  points,  for    v'x  and 

v'y.    are  not  necessarily  proportional  in  those--  zones,  and  dv/dy 

is  equal  to  zero  at  the  axis.    Actual  measurements  of    v'y    (in  air 

with  a  hot-wire  anemometer)  (13)  yield  curve    B    of  figure  9,  which 
is  seen  to  have  a  finite  value  at  the  axis,  to  increase  toward  the 
wall,  and  finally  to  drop  off  abruptly.    The  effect  of  the  unknov/n 
proportionality  factor  is  indicated  by  the  displacement  of  the  two 
curves . 

Nikuradsc  determined  the  trends  of  £  and  £   from  the  slope 
of  measured  velocity-distribution  curves  in  boi-h  smooth  and  rough 
pipes,  as  shovai  by  curves    C    and    G    in  figures  10  and  11.  The 
magnitude  of  Jl   was  assumed  to  attain  a  maximum  value  at  the  axis, 
where  one  would  expect  the  eddies  to  have  the  greatest  size,  al- 
though equation  28  actually  becomes  indeterminate  when    y/r0  =  1 

and    dv/dy  =  0.    The  magnitude  of  £  in  equation  29  becomes  inde- 
terminate for  the  same  reason,  although  £  wil.TL  approach  zero  if 
v'y  does. 

TOiile  equations  3  and  9  were  developed  for  the  velocity  dis- 
tribution in  the  wall  region,  they  have  been  found  to  agree  with 
: measurements  almost  as  far  as  the  pipe  axis.    Differentiation  of 
■either  expression  would  indicate  that    dv/dy    Is  inversely  propor- 
tional to    y  ,  which  yields  curves    D    and    H    of  figures  10  and 

11.  Both  Jt  and  <f  now  reach  zero  at  the  axis  as  well  as  at  the 
wall.  On  the  other  hand,  Von  Xarman's  similarity  hypothesis  for 
turbulence  in  the  intermediate  zone  between  wail  and  axis  led  to 

the  expression  £.  =  k  (dv/dy)/(d  v/dy~),  introduction  of  which 
into  equations  28  and  29  results  in  curves    E    and    I    for  £  and 

£  .    Curve    F    was  derived  in  a  statistical  treatment  of  turbu- 
lence by  Gcbelein  (2). 
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Despite  the  departure  of  these  several  (curves  from  one  an- 
other in  both  form  and  position,  the  corresponding  expressions  for 
velocity  distribution  and  resistance  are  in  sufficient  agreement 
with  measurements  for  a.U  practical  purposes — ceven  v/hen  applied  to 
uniform  flow  in  open  channels  (except  for  possible  adjustment  of 
consta_nts).    On  the  other  hand,  since  the  distribution  of  suspended 
matter  depends  directly  upon  <£  as  a  function  apf    y  ,  one  hesitates 
to  use  any  of  these  distribution  curves  without:,  further  experimen- 
tal check.    Unfortunately,  the  hot-wire  method",  of  measuring  turbu- 
lence in  air  does  not  seem  to  be  suited  to  use  in  water;  and  the 
simulation  of  a  water  surface  with  air  as  the  miodcl  fluid  is  obvi- 
ously out  of  the  question.    Several  alternat j.vtc  methods  are  now  in 
an  early  state  of  development  (4,  S,  3),  but  a,cs  yet  no  usabj.e  open- 
channel  measurements  are  at  hand. ' 

Leighly  (5),  Christiansen  (i),  end  others  have  sought  to 
correlate  the  distribution  of  sediment  and  velocity  in  actual 
streams  by  determining    dv/dy    from  measured  velocity  profiD.es  and 
then  computing  £   from  equation  29.    Such  resmlts,  however,  are  no 
more  dependable  than  Nikuradse's  curve    G  ,  f o:r  the  reasons  already 
stated.    Moreover,  there  is  not  yet  sufficient  evidence  to  show 
whether  the  magnitude  of   <f  determined  from  equation  29  is  numeri- 
cally equal,  or  merely  proportional,  to  that  i"  -  r  the  sediment  sus- 
pension even  in  the  zone  for  which  this  equation  indicates  the 
proper  functional  trend. 

Whatever  the  true  analytic  relationship  'between  £  and    y  , 
certain  general  conclusions  may  be  drawn  by  accepting  one  which  docs 
not  depart  too  seriously  from  probable  conditions.    The  frequent 
assumption  that   <f  is  constant  over  the  entire  flow  section  is,  of 
course,  far  from  the  truth,  although  it  may  be.  come  more  nearly 
constant  some  distance  from  the  boundary  than,  is  indicated  by  the 
curves  of  figure  11.    At  a  free  surface  both  M.  and    v'v    are  un- 
doubtedly  smaller  than  at  the  axis  of  a  pipe,    approaching  zero  if 
the  free  surface  is  smoothj  that  <f  does  not  ilways  approach  zero 
as  a  limit  is  indicated  by  visible  disturbances  at  the  surface  of  a 
very  turbulent  stream. 

Thus,  except  very  close  to  the  bed  and  fee  the  surface,  either 
curve    H    or  curve    I    may  still  indicate  the  [general  trend  of  <f 
in  a  wide  channel.    Selecting    H    as  the  simplcer  of  the  two,  equa- 
tion 2/+  may  be  integrated  to  yield 


2- 


(30) 
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wherein  2  =         -\Z7/p  -  ^l/f/Q     •    In  a  semi-logarithmic 

diagram,  the  relative  position  of  the  distribut'&on  curve  should 
then  be  governed  by  the  magnitude  of  &  ,  and  the  form  by  the  mag- 
nitude of    D/a  ;  the  parameter    a    need  be  only  so  grea.t  as  to 
insure  fully  developed  turbulence  for  all  positive  values  of    y  . 
Figure  12  thus  shows  schematically  a  series  of  relative  distribu- 
tion curves  for  multiple  values  of  2r  and  for    I '/a  -  1  =  10. 

As  has  already  been  remarked,  any  relationship  of  this  form 
will  not  yield  absolute  values  of  sediment  concentration,  for  it 
provides  no  means  of  evaluating  •  c-  .    On  the  etcher  hand,  once  the 

concentration  and  settling  velocity  of  each  graim  size  at  level  a 
are  known,  the  absolute  distribution  curves  for-  each  size  may  be 
obtained.    Thus,  from  a  sample  of  the  suspension  at  the  reference 
level,  one  could  determine  the  settling  velocities  and  the  absolute 
curves  for  each  value  of    w    included  in  the  sasnple,  integrate  the 
product    cv    for  each  curve,  and  then  add  all  together  to  obtain 
the  total  rate  of  transport  of  suspended  load  ai  ovc  the  reference 
level.    This  procedure  would  obviously  be  a  tedious  one,  and  would 
have  to  be  repeated  for  every  change  of  flow  conditions.    One  is 
therefore  inclined  to  seek  a  more  convenient  relationship  in  terms 
of  average  sediment  parameters. 

As  in  the  case  of  the  distribution  of  pariticle  size,  an  anal- 
ysis of  the  fail-velocity  characteristics  of  a  sediment  sample  may 
be  plotted  in  histogram  form  (fig.  13a),  the  abscissa,  scale  repre- 
senting the  logarithm  of    w    and  the  ordinate  scale  the  concentra- 
tion   c    per  class  interval.    The  resulting  block  diagram  will 
approach  a  smooth  weight-frequency  diagram  as  tlae  class  interval 
becomes  small,  the  area  enclosed  by  this  curve  froing  equal  to  the 
concentration    C    of  the  total  ria.terial  in  the  sample.    If  one 
assumes  curve    A    to  represent  the  weight -frequency  distribution  of 
w    at  the  reference  level    a  ,  it  will  be  found',  that  the  frequency 
curves  at  various  levels  above    a    will  be  as  sfhown  in  figure  13b 

for  the  function  <f  =  ft  )  plotted  in  figrarc  12 — the  enclosed 

areas  again  representing  the  magnitudes  of    C    ;at  the  respective 
levels . 

Were    c    and    C    identical  at  every  level,  (that  is,  were  the 
entire  range  of    w    only  one  class  interval),  fit  is  evident  that 
the  vertical  distribution  curves  for    C    would  (coincide  in  figure 

12  with  the  corresponding  curve  for   2  =  w/'c  ynj^p  .    For  a  given 

value  of  k  y'T/p  ,  actual  departures  from  this  limiting  case  will 

depend  primarily  upon  the  geometric  mean  abscissa     wm    of  the 


W- frequency  curve  at  the  level    a  ,  secondarily  upon  the  extent  to 
■which    w    deviates  from  this  mean,  and  to  a  practically  negligible 
degree  upon  the  skevmess  or  asymmetry  of  the  curve.    For  instance, 
curves    A,  B,  and  C    of  figure  13  display  the  same  proportional 
variation  of    w  ,  but  differ  in  geometric  mean;  for  such  frequency 
curves,  at  levels  a__,  _  the._ departure  of ..  the  .respective.,  distribution 
curves  for    C/Ca    from  those  of    c/ca    for  the  corresponding  values 
of    wm    will  be  seen  from  figure  12,  the  difference  evidently /in- 
creasing with  increasing  values  of  y 'a 'y^yp    •    Similarly,  | 

curves    B    and    D  ,  or    C    and    F  ,  have,  respectively,  the  same 
means  but  different  deviations  from  these  means;  the  departure  is 


deviation.    On  the  other  hand,  the  extreme  skowness  of  the  fre- 
quency diagrams    F    and    G    is  seen  to  have  little  effect  upon  the 
distribution  of    C/Ca  ,  the  corresponding  distribution  curves  in 
figure  12  varying  to  a  relatively  small  extent  from  that  for  the 
symmetrical  diagram    C    having  the  some  mean  and  deviation. 

The  magnitudes  of  the  geometric  mean  and"  the  standard  devia- 
tion of  a  frequency  curve  may  be  shown  statistically  to  correspond 
to  the  abscissa  of  the  center  of  gravity  of  the  enclosed  area  and 
to  the  radius  of  gyration  of  the  area  about  a  vertical  axis  through 
this  point.    These  parameters  may  be  taken  directly  from  a  cumula- 
tive plot  of  the  frequency  curves  of  logarithmic  probability  paper. 
On  such  a  plot  of  frequency  diagram,  having  the  symmetrical  form  of 
normal  error  curve,  vd.ll  become  a  straight  lire,  varying  in  posi- 
tion with  the  geometric  mean  and  in  slope  vfith.  the  standard  devia- 
tion.   The  value  of    wra    will  correspond  to  the  intersection  of  the 
cumulative  line  with  the  50  percent  ordinate,  while       ,  the  stand- 
ard geometric  deviation  (i.  c.,  antilogarithm  of  the  standard  devi- 
ation), may  be  found  by  dividing  the  84-.  1  percent  intercept  by  wm 
or    wm    by  the  15. 9  percent  intercept.    For  example,  from  figure  14 
it  is  seen  that  curves    B    and    D  ,  A  ,  and    C    and    E  yield 
values  of    wm    equal  respectively  to  l/8,-/274a  and  1  cm/sec;  like- 
wise, for    A  ,  B  ,  and    C,         =  2,  while  for    D    and    E  ,    <7VI  =  -/2 
In  the  case  of  skewed  curves,  a  close  approximation  may  be  made 
through  use  of  two  terms  of  a  logo.rithmic  Grara-Charlier  series, 
which  yields  a  family  of  curves  coinciding  at  the  Q/^.l  percent  and 
15.9  perceiro  intercepts  with  a  straight  line  having  tho  seme  goo- 
metric  mean  and  standard  deviation.    As  seen  from  curves    F    and  G 
in  figure  14,  the  corresponding  values  of    vjfo    and  may  be  read 

on  the  straight  line  passing  through  these  intercepts. 

From  the  foregoing  discussion  of  the  various  paranetcrs  gov- 
erning the  distribution  of  suspended  load,  it-  follows  that  the 


seen  to  decrease  with  decreasing  values  of 


and  of  the 
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distribution  of    C/C£    over  the  vertical  section  will  depend  upon 
the  three  dimensionless  parameters    a/D  ,  fas,  -|/ T/p      >  an<^  j 

for  a  given  distribution  of    £  .    Since  the  distributions  of  both 
£  and   v    are  governed  only  by  .the  mean  velocity  and  the  relative 
roughness  of  the  bed  in  fully  turbulent  flow,  it  would  appear  that 
the  total  transport  of  suspended  load  might  eventually  be  expressed 
in  terms  of  a  few  pertinent  parameters  characteristic  of  the  bound- 
ary, of  the  flow,  and  of  the  sediment — a  method  of  approach  which 
will  be  discussed  further  in  the  following  section  on  Prediction  of 
Total  Load* 

1 

Before  leaving  the  subject  of  sediment  characteristics,  sev- 
eral, points  will  bear  elaboration.    Despite  the  significance  of  the 
fall  velocity  in  the  study  of  suspended  load,  it  is  general  prac- 
tice to  determine  the  size-frequency,  rather  than  the  w-frequency, 
characteristics  of  sediment  samples.    Sieving  is  used,  so  far  as 
possible,  with  recourse  to  the  pipette  or  other  hydraulic  methods 
for  material  passing  the  200-mcsh  sieve.    It  so  happens  that  sedi- 
ments finer  than  200  mesh  will  generally  settle  in  water  in  close 

accordance  with  Stokes'  law,    w^d    ;  while  the  settling  of  coarser 

materials  will  approach  the  limit    w~d2  as    d    becomes  increas- 
ingly greater  (10).    Moreover,  sieve  analysis  does  not  indicate 
true  hydraulic  diameter  (that  of  a  sphere  having  the  same  density 
and  settling  velocity),  since  it  docs  not  properly  evaluate  shape 
and  density  factors.    Prediction  of  the  w- frequency  from  the  size- 
frequency  curve  is  therefore  at  best  uncertain^  and  because  v-r- 
frequency  and  size-frequency  curves  lying  in  part  beyond  the 
Stokes'  range  will  have  different  degrees  of  skewness,  it  cannot  be 

expected  that    wm    and    <^7    will  always  correspond  to  the  geometric 
mean  and  the  standard  geometric  deviation  of  the  size-frequency 
diagram.    A  direct  method  of  determining  the  v;- frequency  of  coarse 
material  as  well  as  fine  has  still  to  be  perfected. 

In  this  regard,  it  should  be  noted  that  material  finer  than 
a.bout  15  microns  will  almost  Invariably  flocculate,  except  in 
streams  in  arid  climates.    Obviously,  laboratory  analysis  of  de- 
flocculated  samples  may  yield  grading  curves  which  bear  no  relation 
to  the  w-frequency  of  the  actual  suspension.    This  danger  could  be 
obviated  by  analysis  of  suspended  matter  in  the  stream  water  ob- 
tained with  the  sample.    It  might  be  added  that  in  suspensions  of 
appreciable  concentration  the  settling  velocity  of  coarse  particles 
in  the  presence  of  fines  is  less  than  in  clear  water.    The  question 
remains  for  the  present  an  academic  one,  for  the  analysis  of  turbu- 
lence has  not  yet  advanced  to  the  treatment  of  density  variation 
across  the  flow  section. 
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Prediction  of  Total  Load 

Mention  was  made  at  the  end  of  the  section  on  Turbulence  and 
Bed  Load,  of  the  fact  that  the  transport  function  for  bed  load  alone 
might  well  be  extended  empirically  beyond  the  aone  of  appreciable 
saltation  and  thus  used  to  predict  the  total  load  of  a  stream.,  In 
the  same  manner,  one  is  inclined  to  speculate  on  the  possibility  of 
estimating  the  total  stream  load  through  use  of  only  the  function 
for  suspended  load.    The  latter  is  by  far  the  More  venturesome 
method  of  attsck,  but  it  would  seem  at  the  present  time  to  provide 
the  sole  analytical  means  of  expressing  the  transport  capacity  of  a 
stream  in  terms  of  readily  measurable  bed  and  .flow  parameters.. 

Since  such  a  method  would  evidently  require  extrapolation  of 
the  suspended- load  function  well  into  the  bed  '.region,  it  would 
first  be  well  to  note  the  reason  why  this  procedure  has  hitherto 
been  avoided.    If  the  reference  level  in  equation  30  is  placed  at 
the  bed  itself,  a    becomes  equal  to  zero,  and  inspection  of  the 
equation  will  show  that  the  concentration  will  also  be  zero  over 
the  entire  section.     In  other  words,  the  fact  that  £   in  equation 
29  becomes  zero  at  the  wall  would  indicate,  paradoxically  enough, 
that  no  sediment  can  bo  lifted  off  the  bed,  although  an  indefinite 
amount  may  be  held  in  suspension  once  it  rises  above  the  level  of 
zero  mixing. 

"Held  in  suspension"  is  a  phrase  which  bias  a  purely  statisti- 
cal significance,  for  the  process  of  turbulent,  nixing  entails  a 
continuous  interchange  of  material  between  neighboring  levels  of 
the  flowc    The  presence  of  material  in  suspension  therefore  indi- 
cates not  only  that  it  was  originally  lifted  .from  the  bed  in  one 
way  or  another,  but  that  some  means  continues,  to  exist  whereby  new 
material  may  be  carried  into  suspension  to  replace  that  constantly 
settling  out.    In  other  words,  either  t    does  not  reach  a  value  of 
zero  at  the  bed,  or  else  the  momentary  fluctuations  from  a  statis- 
tical mean  value  of  zero  arc  sufficient  to  provide  the  initial 
upward  impulse. 

With  reference  to  figure  3c,  it  is  apparent  that  in  the  ab- 
sence of  viscous  influence  the  bed  roughness  sidy  be  considered  to 
produce  local  disturbances  depending  in  size  "upon  the  parameter  k 
and  in  intensity  upon  the  mean  velocity  at  the  distance    k  from 
the  mean  bed  level.    While  these  eddies  may  actually  yield,  a  zero 
magnitude  of   6  at  zero  elevation  if  evaluated  in  a  statistical 
sense,  the  local  effectiveness  of  the  nixing  process  surely  extends 
even  belo\7  the  mean  level  of  the  bed.    It  doers  not  seem  unreason- 
able, therefore,  to  presume  that  £  remains  es-sentially  constant 
from    y  =  k    to    y  =  0. 
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Under  these  circumstances,  k    would  servre  as  the  Uniting 
value  of    a    in  equation  30.    The  magnitude  of  night  then 

be  approximated  from  equation  26,  using  for  ihve  constant  magnitude 
of  £  that  at    y  =  k  ,  and  for    C0    the  concentration  of  the  bed 
material.    It  is  barely  possible,  moreover,  thxat    a    might  be  set 
equal  to  some  constant  fraction  of    k  ,  Ca    thten  being  taken  equal 
to    C0    and  the  intermediate  step  using  equati'.on  26  being  omiitcd 
completely. 

The  latter  assumption  is  obviously  a  bold  one.    In  fact, 
either  procedure  must  be  subjected  to. thorough!  experimental  check 
before  it  can  be  considered  valid..    Experimental  measurements  (would 
be  necessary  in  any  case,  moreover,  to  determine  the  proper  numeri- 
cal value  of    k  .    Nikuradse,  it  will  be  recaU.cd,  adopted  the 
grain  diameter  as  a  measure  of  roughness,  althiough  the  effective 
height  of  the  irregularities  in  the  sand  coatahng  that  was  applied 
to  the  pipe  was  not  necessarily  the  same  as  tlhe  diameter  of  the 
grains  themselves.    Shields  also  used  the  grcdin  diameter  as  a  meas- 
ure of  roughness  in  the  case  of  initial  movement ;  as  bed  undula- 
tions develop  in  size,  her/ever,  the  individual.,  sediment  grains  no 
longer  determine  the  absolute  bed  roughness^  li:    then  varying  with 
the  size  and  shape  of  the  undulations. 

Equation  12}  with  possible  modification  of  constants,  might 
still  be  used  to  estimate  the  effective  magnitude  of    k    for  an 
open  channel,  for  this  relationship  has  both  analytical  and  experi- 
mental foundation.    Engineers,  on  the  other  ha\nd,  are  more  accus- 
tomed to  thinking  in  terms  of  the  Chezy  equation 


(3D 


in  which    C    is  inversely  proportional  to  the  'square  root  of  the 
resistance  coefficient,  f  : 

C  -l/M  (32) 
f 

C    may  therefore  be  evaluated  either  in  terms-:  of  the  relative 
roughness  of  equo.tion  12,  or,  as  is  more  customary,  through  use  of 
the  empirical  formula  of  Manning: 

/  ^/?a'    /P  6 

(33) 

Although  it  is  not  generally  realized,  the  Manning  n  is  the  sixth 
root  of  an  absolute  roughness  parameter,  n  representing  a 


18 


relative  roughness  similar  to  the    d/r0    used  by  Nikuradse.  The 
Manning  relationship  therefore  applies  to  flovr  In  which  the  in- 
fluence of  viscosity  is  negligible.    If  equation  33  were  truly 
valid,  n°    would  be  directly  proportional  to  values  of    k  =  d  de- 
termined from  equation  12.    The  actual  relationship  may  be  seen 
from  figure  15,  which  may  be  used  to  estimate  at  least  the  order 
of  the  magnitude  of  the  absolute  roughness  from  measurable  flow 
parameters . 

I 

Whether  or  not  the  foregoing  means  of  predicting  total  load 
proves  feasible,  the  method  of  analysis  suggests  a  systematic 
organization  of  the  variables  belonging  in  a  general  transport 
function o    It  was  shown  in  the  section  on  Turbulence  and  Suspended 
Load  that  the  distribution  curve  for  suspended  load  is' a  function 
of  the  material  characteristics  at  the  reference  level,  the  height 
of  the  reference  level  in  its  ratio  to  total  depth,  the  friction 
velocity,  and  the  distribution  of   £  .    Setting    a    equal  to    k  , 
these  parameters  may  be  grouped  dimonsionlessly  in  the  functional 
form 


(34) 


Quantity  in  suspension  =  ^ }         3  J/J^Z^  f 

the  suspension  increasing  with  the  magnitude  of  each  dimensionless 
term.    The  form  of  the  function  depends,  evidently,  upon  the  dis- 
tribution of  <f   ;  but,  once  established,  the  variation  of  E  AZ^^/jyp 
should  remain  a  unique  function  of    y/D  alone.. 

The  magnitude  of         ,  in  turn,  should  vary  with  the  bed 
concentration,  the  absolute  roughness,  the  mean  fall  velocity  and 
its  standard  geometric  deviation,  and  the  effective  nixing  at  the 
bed,  that  is, 


(        /r  /i 


-    .  cr  }  (35) 


If  £  is  evaluated  in  terms  of  k  and  v^-  ,  v^  being  found 
from  equation  9,  it  will  then  be  seen  that 


Q  -  KfC  ,  HZP  ,  <r  (36) 


of    C\:    is  used  in  equation  34a-  ^nG  form  of  the 
function  will  change,  but  except  for    C0    the;  parameters  will  be 
the  same  as  before. 

For  fully  effective  boundary  roughness  :Lt  has  been  shown 
that  the  velocity  distribution  is  a  function  ©if  the  friction 
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velocity  and  the  relative  roughness,  both  parameters  already  ap- 
pearing in  the  foregoing  relationships.    It  sfchould  therefore  be 
possible  to  express  the  total  material  transported  in  the  form 


L      V  T/P       cr  )  (37) 


J  J  tv 


This  expression  for  total  load,  evidently,  iss  the  counterpart  of 
the  empirical  bed-load  equation,  likewise  requiring  experimental 
determination  of  the  functional  form  and  the  .'magnitude  of  the 
numerical  constants.    However,  the  resulting,  ^relationship  will  not 
be  entirely  empirical  because  of  the  choice  amd  arrangement  of  the 
variables  by  analytical  means.    If  it  is  trul;:y  general,  on  the 
other  hand,  it  should  describe  the  earlier  stages  of  movement  at 
least  as  well  as  does  the  bed-load  function.     As  a  matter  of  fact, 
when  simplified  to  the  form 


(33) 


-77 

it  will  be  found  to  be  essentially  the  same  aas  the  Shields  expres- 
sion (equation  20),  since    D  —  l/S  ,  and    Yf/k-.~(/f         )d  within 
the  limited  range  of  size  which  Shields  studied. 

Like  the  bed-load  equation,  the  general?,  relationship  can  be 
rewritten  in  other  forms  through  introduction!  of  the  Chczy  expres- 
sion and  the  equilibrium  equation  for  boundary  shear,    7"  -  fD£>. 


In  particular  it  should  be  noted  that  y  7~/p     -  V 'j/^/s  a  while  f 

in  turn  is  a  function  of  k/D  •  thus  the  fri'dtion  velocity,  y^/ps 
may  be  replaced  by  the  mean  velocity  of  flow..  V  ,  elimination  of 

f    merely  changing  the  form  of  the  general  fumction.    Equation  37 
might  then  be  written  in  the  form 

Although  this  general  function  for  totr/i.1  load  may  be  valid 
even  for  relatively  high  sediment  concentration,  it  must  be  noted 
that  under  such  conditions  the  boundary  roughincss  may  not  be  fully 
effective.    That  is,  the  mechanical  friction,  .resulting  from  con- 
tact of  the  sediment  grains  in  the  zone  of  high  concentration  near 
the  bed  will  tend  to  reduce  the  turbulence  of"  the  boundary  flow — 
producing  in  the  extreme  case  a  region  of  n on-:- turbulent  movement 
at  the  bod  rcughly  equivalent  to  the  laminar  'boundary  layer — 
mechanical  friction  now  taking  the  place  of  f?luid  viscosity.  This 
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circumstance,  coupled  with  the  failure  of  existing  equations  for 
velocity  and  sediment  distribution  to  take  into  account  an  appreci- 
able density  gradient,  makes  the  problem  of  high  sediment  concen- 
tration decidedly  more  complex  than  the  conditions  herein  assumed. 

Conclusion 


Analysis  of  sediment  transportation  in  the  light  of  fluid 
turbulence  at  once  emphasizes  the  need  for  experimental  data,  both 
to  provide  the  necessary  numerical  constants  for  proved  functional 
relationships  and  to  test  the  validity  of  still  unproved  assump- 
tions.   One  reason  for  the  lack  of  such  data  has  lain  in  the  dif- 
ficulty of  systematizing  research  without  some  immediate  experimen- 
tal goal'.    It  is  to  this  end  that  the  present  analysis  has  been 
prepared,,  for  it  provides  a  definite  program  of  research  in  the 
proof  or  disproof  of  methods  herein  proposed.    By  way  of  summary, 
therefore,  it  "would  seen  expedient  to  outline  the  points  in  this 
analysis  still  awaiting  experimental  study,  with  suggestions  as  to 
suitable  procedure. 

(1)  Shields'  treatment  of  initial  movement  and  subsequent 
transportation  of  material  of  uniforn  size  and  shape  merits  early 
investigation  as  to  the  effect  of  change  in  size  and  grading. 
Logical  parameters  for  the  arbitrary  control  of  size  distribution 
are  the  geometric  mean  diameter  and  the  standard  geometric  devia- 
tion of  size-frequency  curves.    Shape  is  probably  of  secondary 
importance  in  natural  bed  material. 

(2)  The  problem  of  open-channel  flow  is  replete  with  unknown 
factors,  even  for  channels  of  very  wide  cross  section.    The  loga- 
rithmic velocity-distribution  and  resistance  relationships  must  be 
checked  systematically,  and  numerical  constants  determined  for 
cases  of  actual  as  well  as  artificial  roughness.    The  distribution 
of    v'y    ,  j£    ,  and  particularly  <f  must  be  determined — with  ap- 
paratus still  to  be  perfected. 

(3)  Even  prior  to  the  direct  measurement  of  £  as  a  function 
of    y  ,  the  distribution  of  suspended  load  may  be  correlated  with 

"■  the  velocity  distribution  through  controlled  laboratory  tests.  It 
would  be  most  expedient  to  avoid  the  dependence  of  the  suspension 
upon  bed-sediment  roughness  by  securing  artificial  roughness  ele- 
ments (metal  projections)  to  the  flume  bottom  in  such  a  manner  that 
the  bed  material  could  circulate  freely  without  altering  the  magni- 
tude of  the  relative  roughness.    While  this  investigation  might 
best  begin  with  low  concentrations  of  uniform  material,  the  nature 
of  the  artificial  roughness  is  such  as  to  permit  studies  of  the 
effects  of  grading  and  of  high  concentration  upon  the  suspension  in 
the  simplest  manner  possible. 
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(4)    Dcterr.iinc.tion  of  stream  capacity  may-  proceed  in  two 
ways.:    First,  the  relationship  nay  be  sought  beStween  the  suspended 
load  and  conditions  at  the  bed,  by  investigating  the  possibility  of 
evaluating  the  sedinent-distribution  .function  iln  terns  of  the  bed 
roughness  end  the  concentration  of  the  bed  naferial.    Secondly,  and 
of  greater  practical  importance,  a  general  funcrtion  of  the  usual 
bed-load  type  nay  be  extended  to  include  total  load  by  continuing 
experiments  beyond  the  stage  of  appreciable  sanitation;  evidently, 
the  usual  bed-load  trap  will  not  catch  materia!!!  in  suspension,'  a 
vertical  sampler  of  the  divisor  type  being  a  possible  substitute; 
whatever  the  experimental  technique,  the  investigation  must  deter- 
mine the-  functional  relationship  of  dimens ionises s  parameters 
similar  to  those  in  equation  39. 
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Tabulation  of  Symbols 


a  =  reference  level  for  suspended-load  function. 

c  =  sediment  concentration  by  weight  for  given  fall  velocity. 

d  =  dianeter  of  sedinent  grains . 

e  =  base  of  natural  logarithms. 

f  =  coefficient  of  resistance  to  flow. 

g  =  acceleration  of  gravity. 

k  =  absolute  roughness  of  boundary. 

1  =  length  scale  of  turbulent  nixing. 

n  =  Manning  roughness  coefficient. 

r0  =  pipe  radius. 

v  =  mean  velocity  of  flow  at  any  point » 

V'x  =  ■  velocity  of  fluctuation  in  direction  of  flow. 

v'y  =  velocity  of  fluctuation  normal  to  boundary. 

w  =  settling  velocity  of  sedinent. 

Vfm  =  geometric  nean  settling  velocity. 

y  =  normal  distance  from  boundary. 

z  -  exponent  in  the  suspended-load  function. 

C  =  total  concentration  by  weight  at  any  level  (also  used 
locally  for  Chezy  coefficient) . 

D  =  total  depth  of  flow. 

F  =  force  exerted  by  flow  on  bed  material. 

G  =  rate  of  sediment  discharge  by  weight. 

Q  =  rate  of  total  discharge  by  weight. 
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R      =  hydraulic  radius  (also  used  locally  for  resistance  of  bed 
material  to  initial  novenent) . 

f?     =  Reynolds  nur.iber. 

S     =  slope. 

T      =  tractive  force  (intensity  of  shear  at  boundary). 

Tc    =  tractive  force  at  beginning  of  bed  novenent . 

V     =  mean  velocity  of  flow  for  cross  section. 

u    =  coefficient  for  sediment  characteristics . 

p   =  proportionality  factor  for  nixing  process. 

<J*  =  specific  weight  of  fluid. 

J*  =  specific  weight  of  -sedinent . 

8  =  thickness  of  boundary  layer. 

€  =  coefficient  of  turbulent  nixing. 

■5  -  coefficient  of  sedinent  resistance  to  flow. 

ft  ~-  flow  constant,  equal  to  about  0.4-0. 

M  =  dynanic  viscosity. 

V  -  kinenatic  viscosity. 

P  -  nass  density. 

^  =  standard  geonctric  deviation  of  \?-frequency . 

I  -  intensity  of  fluid  shear . 

2~  -  intensity  of  shear  at  boundary. 

P  =  function  of. 
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